Boundaries of Subcritical Coulomb Impurity Region in Gapped Graphene 
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The electronic energy spectrum of graphene electron subjected to a homogeneous magnetic field 
in the presence of a charged Coulomb impurity is studied analytically within two-dimensional Dirac- 
Weyl picture by using variational approach. The variational scheme we used is just based on utilizing 
the exact eigenstates of two-dimensional Dirac fermion in the presence of a uniform magnetic field 
as a basis for determining analytical energy eigenvalues in the presence of an attractive/repulsive 
^~ > ' charged Coulomb impurity. This approach allows us to determine under which conditions bound 

state solutions can or can not exist in gapped graphene in the presence of magnetic field. In addition, 
the effects of uniform magnetic field on the boundaries of subcritical Coulomb impurity region in 
the massless limit are also analyzed. Our analytical results show that the critical impurity strength 
decreases with increasing gap/mass parameter, and also that it increases with increasing magnetic 
field strength. In the massless limit, we investigate that the critical Coulomb coupling strength is 
independent of magnetic field, and its upper value for the ground-state energy is 0.752. 
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Since the discovery of graphenei , strictly two-dimensional (2D) carbon system with hexagonal lattice structure, a 
great deal of both experimental and theoretical research efforts has been achieved to identify the electronic structures 
of graphene and graphene based nanostructures such as graphene dots and graphene nanoribbons^. In fact, from the 
theoretical point of view, after Wallace's work 3 , Semenoff 4 was the first to construct the relativistic 2+1 dimensional 
electrodynamics analogue of tight-binding model of electrons in a 2D hexagonal lattice. In other words, he investigated 
that, in the long- wave length (continuum) limit, the low-energy electronic structure of this model exhibits linearly 
I ■ dispersing massless fermionic quasi particles around the degeneracy points obeying to 2+1 Dirac- Weyl Hamiltonian. 
As a consequence, all the single particle low-energy physics of graphene is believed to be well simulated by the 
Dirac- Weyl equation. 

Due to their relevance for the transport properties, impurity effects are, at present, at the center of both theoretical 
and experimental investigations in graphene physics^—. There are non-negligible effects especially on the electronic 
properties of graphene based devices, since the Dirac points of undoped graphene which are very sensitive to such 
effects^ 3 -. Additionally, a recent experiment 1 ^ performed by depositing potassium atoms onto undoped graphene 
addresses charged Coulomb like behavior. Moreover, it is shown that a pronounced asymmetry occurs in transport 
cross-section depending upon Z — > — Z for the Coulomb scattering^— , and it is also believed that the long-range 
t-H ■ charged Coulomb impurity limit the mobility of graphenei 14 ' 19 . 

From the quantum electrodynamical (QED) point of view, it is well-known that the interactions between the charged 
I . particles are associated by the exchange of virtual photons. These interactions are represented by a fermionic Coulomb 
^-H ' potential which is identified with its perturbative expansion in powers of the fine structure a = e 2 /he — 1/137. 
0^ . However, it is broken down in the case of superheavy nuclei, atoms or quasi-molecules where a is of order of unity 2 ^. 
1 This hypothetical regime with Za > 1 which requires a non-perturbative treatment (due to difficulties of convergence 
of expansions) has been of academic interest in the past. However, due to the fact that the effective coupling constant 
of graphene is at order of a = e 2 jehvF ~ 1 for conventional Si02 substrate^— ' 15 ' 16 , condensed matter analog of this 
regime is of experimental interest, and it is now at the center of theoretical investigations. In particular, "mass" or 
energy gap in graphene is frequently interpreted as being related to chiral symmetry breaking of 2D massless Dirac- 
Weyl fermions induced by substrate 2 ^, spin-orbit coupling 22 , and by boundary conditions^ 3 -. Since this enables one 
to tune a gap by just adjusting the external parameters of graphene, it is of particular interest itself. 

In the present paper, we investigate analytically the effects of a Coulomb impurity with effective charge Ze onto 
the energy spectrum of 2D massive fermions, i.e., gapped graphene, with a uniform magnetic field perpendicular to 
the graphene plane. In other words, we analyze the effects of uniform magnetic field onto the graphene electron 
with impurity, and hence we determine the boundaries of subcritical regime. As mentioned above, due to both 
strong coupling and long-range characters of the Coulomb impurities in graphene, these effects can not be addressed 
by conventional perturbative techniques. Therefore, we suggest a variational approach to account for such effects, 
induced by the impurity. 

The Dirac Hamiltonian with electromagnetic potentials that we describe the motion of an electron and a hole in 
2D graphene subjected to a constant uniform magnetic field, perpendicular to graphene plane, can be written in the 
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H = hv F a ■ (—idp+^r-A^ + f3mv F — eA (r), 



(1) 



where Ao and are the time and space components of the four- vector potential (^4oj A^) , due to impurity Ao = Ze/er 
and static uniform magnetic field with A = B(—y,x) /2, respectively. Here, while Z > refers to the attractive 
impurity potential that binds electrons to impurity and repels holes, or vice versa if Z < 0. In Eq. ([1]), we have 
used the Dirac-Pauli representation of Dirac matrices, ot and (3, each of which are written in two by two-block form. 



Using two component spinor representation as W — {(p" X~ ), 
H^f = E^> satisfies the following coupled first order equation: 



we see that each component of the eigenvalue equation 
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where we write energy, coupling strength and "mass" terms in units of Hvf , as E = E/Hvf, Ao(r) = Ze/hvpr, and 
Mq = mvp/h, respectively. From our knowledge on planar relativistic and non-relativistic electron systems in the 
presence of both Coulomb and magnetic field, it is not possible to find exact analytical solution s 24 ' 25 of Eq. © due 
to having hidden sii algebraic structure^. 

Decoupling Eq. ([2]) in the absence of impurity potential one can easily see that the upper component of 'J, i.e., 4>, 
should satisfy the second order equation 
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whose solutions can then be written in terms of Laguerre polynomials 27 
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provided that the energy quantization condition 



[Av + 2 (|m| + 1) + 2 (s + m)] ^~ = E 2 - M 2 , (5) 

_l_ 1/2 

is fulfilled. It is obvious that Eq. ^ yields infinitely degenerate energy eigenvalues E n = ± [Mq + (2neB / he)] 
with n = v+ [(|m| + m + s + 1) /2] . Inserting Eq. Q into Eq. (|2|), the other component of the spinor can explicitly be 
determined. Therefore, to the positive energy spinors in the absence of impurity potential, we can write the complete 
and orthonormalized solutions of Eq. (fTJ) for s = +1 and s = — 1 cases as 
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and +1 (jp) = _ x fap) — l)"^ (^ 2 p 2 ), where Ll" 1 ' are the well-known associated Laguerre polynomials . In 
Eq. ([6]) and Eq. (JTJ), we have also defined 



N+ m (M ,B) = 



v\ Mo + [M 2 + (2n/£ 2 )] 
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and 



(v + \m\)\ 2[M 2 + (2n/£ 2 )} 1/2 
(M , n) = ^= {Mo + [Ml + (2n/£ 2 )] 1/2 } , 



where -f 2 = eB/2hc, and £ is the magnetic confinement length given by £ = y/hc/eB. Analogously, one can follow the 
same treatment for the negative energy spinors. 

At this work, we suggest the states given by Eq. ^ and Eq. ([7]) as the trial states for the whole system, i.e., 2D 
Dirac equation with charged Coulomb impurity given by Eq. ([2]). Therefore, for the expectation value of H given by 
Eq. ([IJ, we obtain 
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where Ao(r) = Za/r, and a is given by a = e 2 /ehvF- By using first Eq. ([6]) together with Eq. (|7|) in Eq. ([8]), and 

then performing the necessary integrals^, the total variational energy E nms (7) of the whole system can finally be 
written in the form 
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Minimizing Eq. ([9]) with respect to variational parameter 7, we obtain 
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where we restrict ourselves to the physical positive sign of the variational parameter, due to the sign convention we 
employed in the normalizatio n 2 ^ 30 . Here, it should also be noted that, in the literature, different sign conventions 
may be choosen, however some of them may introduce artificial difficulties, especially, in the massless limit ( Ref. 
for a detailed discussion). Therefore, by labeling the electron and hole states A = +1 and A = — 1, respectively, we 
obtain an analytical result for the whole energy spectrum of the system with Coulomb impurity potential as 
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with 
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for all values of Z and Mq if and only if H (Z, Mq) < 1, which implies the condition 
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should be satisfied in order that there might exist oscillator like discrete energy spectrum. Otherwise, Eq. ([12]) becomes 

imaginary. In U(Z,M ) = limit, Eq. (HJ) becomes E n = ±^ M$ + 2n/£ 2 . Thus, we recover the usual Landau 
picture of 2D Dirac equation with mass. 

It should be noted that E n given by Eq. (TT21 depends strongly on the sign of Z, i.e., it is asymmetric with 
respect to reverse of Z, due to the fact that, in the presence of a uniform magnetic field, electron (hole) in a attractive 
impurity potential (Z > 0, [Z < 0)) does not have the same energies as electrons (holes) in a repulsive impurity 
potential (Z < 0, (Z > 0)). This is apparent from FIG. [1] (a) and (b) where the evolution of the ground state energy, 

Ei of massless graphene electron(hole) in the magnetic field is shown for (a) various values of Z > 0, and (b) for 

various values of Z < 0. As is seen from the figure, when Z approaches to ±Z cr , E^ is dramatically reduced, and 
it vanishes at Z = ±Z cr . It can be easily seen from Eq. (p~3|) that, in the case of massless graphene, i.e., Mo = 0, this 
critical value of Z, i.e., Z cr is given by 

Z cr (0,i?)=±^=-— — (14) 

which is independent of i , i.e., of B, and it leads to Z cr — = 0.752 for the ground-state of the system, i.e., for 

n = 1 corresponding to {v, m, s) — (0, 1, —1), and (0, 0, +1). In fact, it is well-known that, in the absence of magnetic 
field, this critical value is equal to 0.5^i& for the massive case. This difference arises from the use of oscillator states 
as trial states in our problem to obtain the dependence of critical value of Z on magnetic field. Therefore, it is obvious 
that our analytical results are valid for high magnetic field regime. 

In the presence of a gap/" mass" term, as shown in FIG. the picture is quite different depending on the values 
of Z, compared with those of found in FIG. [1] (a) and (b). In this figure, we use a gap of value ~ 0.26 eV which is 
reported in single layer graphene due to SiC substrate by Zhou et al^i, and we choose M = O.lt which corresponds 
to 270 meV with hopping energy t — 2.7 eV. Curves for 
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without magnetic field but with impurity^, and curves for E n (Mq,B) = =f^/Mq + (2n/£ 2 ) with magnetic field but 
without impurity are also plotted in the same figure for comparison purposes. Here, sgn is the sign function. We 
represent the former group, i.e., electron and hole energy bands, by horizontal solid lines, while the bold straight lines 
are used for the later one. We first notice that, with inclusion of gap/" mass" term, the behavior of the ground-state 
energy changes from a square root of B dependence to a linear dependence on B. As B is decreased, the curves 
for electron energies with Z = —0.1 and —0.2, and the curves for the hole energies with Z = 0.1 and 0.2, all they 
approach E n [Mq, Z) indicated by horizontal lines which correspond to impurity energy in the absence of magnetic 
field. It should also be noted that, while the curve for gapped-graphene electron (hole) energy with Z = 0.05 (—0.05) 
spreads very little through the gap, the curve for gapped-graphene electron (hole) energy with Z = 0.2 (—0.1) lies 
significantly below the gap. 

To understand this picture better, it is necessary to study the impurity binding energy, which is a measure of 
how much of E n consists of the charged Coulomb impurity. In other words, it is defined as the energy difference 
between the energy of graphene (massless or massive) without impurity and the energy of graphene with impurity, 
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i.e., AE n — ±\/Mq + (2n/£ 2 )— E n . In FIG. [3](a-c), impurity binding energies for low- lying gapless and gapped- 
graphene states are given for three different values of magnetic field strength. From these figures, we see that 
enhancement in magnetic field strength leads to more binding. Moreover, in the absence of gap term, binding energies 
for electrons (thin lines) and holes (thin dashed lines) have the same magnitude, i.e., they are all equal in magnitude 
under Z — > — Z. 

To investigate the influence of gap/" mass" term in detail, we have also examined the dependence of Z cr on both 
magnetic field strength B and gap/mass term Mo, in FIG. BJa) and (b), respectively. The curves of the figure 

corresponds to e\ case only, since tendencies of changing of E x with B and Mq are similar to those of E 1 . For a 
given value of Mq, as B is increased, the minimum value of Z cr curves shifts to larger values. As expected, all figures 
in FIG.|4ja) and (b) illustrate that, as Mq approaches to zero, threshold value of Z cr is given by 0.752 as is calculated 
above. By switching a gap, Z cr is drastically reduced compared to that found in the massless limit. 

In FIG. [SJa) and (b) for comparison purposes, we plotted the evolution of Z cr as a function of (ym) for massless 
and massive cases, respectively. The energy levels are dependent on the principal quantum number n, and they are 
degenerate with respect to to, —to and s. By comparing the left and right panels of FIG. Ufa), we see that the 
degeneracy of energy levels with respect to m, — to are partially removed by the impurity. However, inclusion of a 
gap/"mass" term ( FIG. Efb)) splits the degeneracy of these levels with respect to to, —to and s. Additionally, by 
switching a gap/" mass" term, we see that critical values of Z are (i) smaller than those presented for the massless 
limit, and (ii) they diminish when increases for (v, to < 0,+) while Z cr for (v, m < 0,—) states increases. These are 
due to the fact that, while the massless Dirac-Weyl Hamiltonian with external fields commutes with helicity operator, 
i.e., it is invariant under the chiral transformation, the mass term is not, therefore it breaks the chiral symmetry. For 
the physical origin of this axial anomaly, we refer to the book of Huang— where a detailed discussion of nonexistence 
of conserved axial vector current is given in a whole chapter, and we also refer to the paper of Lee et al^ wherein 
they concluded that first chiral states does not decouple, even in the massless limit, and second massless spinor 
electrodynamics is a pathological theory and therefore it should always be considered to be the limit of a massive 
theory. 

In conclusion, we have investigated the effect of charged Coulomb impurity on the Landau level spectrum of 
both gapless and gapped-graphene near the Fermi point, by using their 2D continuum massless and massive Dirac 
descriptions, respectively, in the framework of a variational procedure. Since, we know that the success of variational 
treatment strongly depends on the choice of trial wave functions together with the set of parameters included, we 
suggest the basis sets of unperturbed Hamiltonian as trial wave functions with one parameter so as to reproduce the 
exact analytical results in the absence of perturbation. Indeed, to test the results of our variational calculation, it is 
enough to look at the change of variational parameter as a function of impurity strength Z. In the limit Z — ¥ 0, the 
variational parameter becomes 7 = l/£. This is an indication of the adequacy of the trial wave function, since it yields 
the exact results in the absence of impurity, as expected. Furthermore, in the absence of impurity the variational 
parameter 7 = \jt is the inverse order of magnitude of magnetic confinement length. When impurity is switched 
on, renormalization of the magnetic confinement length should be taken into account so as to adapt itself to the 
response of the presence of impurity, for which the energy develops a new minimum. This is achieved by adjusting the 
variational envelope wave function with a variational parameter 7 until the energy is minimized. Thus, the inverse 
of 7 may be interpreted as the amount of effective displacement of the magnetic confinement length caused by the 
impurity. Although the results of variational calculation we use give good agreement with the exact ones in the 
absence of impurity, and they also yield exact analytical results for the ground- and exited-states in the presence of 
impurity, such a description may not be accurate for the Landau bands with high indices in the case of high magnetic 
fields. In this case, adjacent Landau bands may overlap. As a consequence, a more sophisticated trial wave function 
composed of linear superpositions of trial wave functions we use with more than one variational parameter may be 
needed for a numerical approach^ 

In summary, by employing a variational procedure based on a choice of trial wave functions as the basis sets of 
Landau levels, we obtained analytical results for both energy eigenvalues and the critical Coulomb coupling strength 
as functions of both magnetic field strength and gap parameter for a graphene electron subjected to a homogeneous 
magnetic field in the presence of a charged Coulomb impurity. The analytical results that we obtained here show 
that (i) the critical impurity strength is independent of the magnetic field strength in the massless limit, and (ii) it 
drastically reduces when a gap/" mass" term is switched on. However, enhancement in magnetic field strength leads 
to increase in critical Coulomb coupling strength. 
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FIG. 3: (Color online) Impurity binding energies of gapless graphene electron (thin solid lines) and gap less graphene hole (thin 
dashed lines) as a function of Z for the three low-lying Landau levels, i.e., (a) for (10), (b) for (20), (c) for (21), for three values 
of magnetic field, B = 1 (black), 3 (red) and 5T (blue). Binding energies of gapped graphene electron and hole states are also 
shown in the same graphs but by thick counterparts for Mo = Mo/t = 0.1. 
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FIG. 4: (a)The critical value of Z cr given by Eq. (|13p as a function of magnetic field for four different gap values, and (b) as a 
function of gap parameter, but for four different magnetic field values. 
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FIG. 5: (Color online) Dependence of the critical value of Coulomb coupling strength, i.e., Z cr on quantum numbers (ym) for 
(a) without gap, A/o = (b) with gap, Ma — 0, 1, respectively. The dashed lines are drawn to guide the eye. In these figures, 
left panels corresponds to m > 0, while the right ones refer to m < 0. 



